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On) , ABSTRACT 

I I ' Context. Determining the parameters of massive stars is crucial to understand many processes in galaxies and the Universe, since 

' these objects are important sources of ionization, chemical enrichment and momentum. 10m class telescopes enable us to perform 

, detailed quantitative spectroscopic analyses of massive stars in other galaxies, sampling areas of different metallicity. Relating the 

• ■ stars to their environment is crucial to understand the physical processes ruling their formation and evolution. 

' Aims. In preparation for the new instrumentation coming up on the GTC, our goal is to build a list of massive star candidates in 

Mh' the metal-poor irregular galaxy IC 1613. The catalogue must have very high astrometric accuracy, apt for the current generation of 

Q , multi-object spectrographs. A census of OB associations in this galaxy is also needed, to provide important additional information 

5_j . about age and environment of the candidate OB stars. 

' Methods. From INT-WFC observations, we have built an astrometric and photometric catalogue of stars in IC 1613. Candidate blue 

' massive stars are preselected from their colors. A friends-of-friends algorithm is developed to find their clustering in the galaxy. While 

' ' , a common physical origin for all the members of the associations cannot be ensured, this is a necessary first step to place candidate 
OB stars in a population context. 

I ' Results. We have produced a deep catalogue of targets in IC 1613 that covers a large field of view. To achieve high astrometric 

^ ' accuracy a new astrometric procedure is developed for the INT-WFC data. We have also built a catalogue of OB associations in 

lO , IC 1613. We have found that they concentrate in the central regions, specially in the Hii bubbles. The study of extinction confirms 

lO ■ that it is patchy, with local values of color-excess above the foreground value, 

'sj" ■ Conclusions. 

' Key words. Stars: early-type - Galaxies: individual: IC I6I3 - Galaxies: photometry - Galaxies: stellar content - Astrometry - 

' Catalogues 
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1. INTRODUCTION figurations since exposure times are long. Targets are often se- 

, , . _ , r , ■ , • , lected based on their magnitude (favouring the brightest can- 

Massive stars are powerful agents of galactic evolution because ^.^^^^^^ ^^^^^^^ ^^^^ ^^^^^.^^^ constrained by the multi- 

of their stellar winds and UV radiation. However, these two object spectrograph used for the ob servations (see for instance 

^ properties depend strongly on the chemical composition of the [Bresolin et al.ll2007l) : ICastro et all (l2008l) ). However, more in- 
. - object, making studies of massive stars in Loca Group galaxies ^^^^idates can be found if the target selection criteria 

of non-solar metallicity compelling. Very metal-poor irregulars ^^^^^ photometry or location but also on the nur- 

have an interesting add-on, local cosmology. By studying their j^tion. Moreover, OB associations are an ideal test 

massive stars and OB associations we can learti how star forma- ^^^^^j^ ^^^^ evolution of massive stars because of their 
tion and evolution works at low metallicity, and infer a template ^mall age spread and homogeneous chemical compo- 

for the deeper Universe. sition (see e.g. Simon-Diaz 2009, in prep.). 
To exploit the improved capabilities that will be available 

with the Gran Telescopio Canarias (GTC), we are building a t .u- i . j n • • i t/-' i/cn /t 

' . ^. , . r , , , In this work we study the magellamc UTegular IC 1613 (Im 

list of candidate blue massive stars in galaxies of the Local t,, , r^^-r — p rp i, j i , r,iA >\ u * • /j- » 

^ ^ . , . ^ ° . . , , BV, IFisher & TuUvl ( 11975ft ). because of its proximity (distance 

Group Spectroscopic observations of massive stars outside the ^^^^^j^^ DM=24.27 lOolphin et al.1 JMl ). low stellar density 

Milky Way and the Magellanic Clouds require multi-obiect con- j u- ui i j iT i * j • 

J J ^ 1 J and highly res olved pop i ilation , very low fore ground reddening 

Send offprint requests w.M. Garcia (e-mail: mgg@iac . es) (E(B-V)=0.02EiiitiD CHI, 0.03iiidia|r971|) ) and most 

Based on observations made with the Isaac Newton Telescope op- importantly its low metallicity. The Study of the emissionjme 



Dciscu uii uusci vaiiuiis iiinuc Willi Liic isatiu iNcwiuii icicsuupc up- " - - I 

erated on the island of La Palma by the Isaac Newton Group in the SEgCtra of some Hll regions yielded 0.04 to 0.08 dTalent, 



Spanish Observatorio del Roque de los Muchachos of the Instituto de 1980; Davidson & Kinman, 1982; Dodorico & Dopita, |1983 

Astrofisica de Canarias. Peimbert et al., 1988; Kingsbi 

** Based on observations made with ESO Telescopes at Paranal lar abun dances of B supergia n 

Observatories under programme ID 078.D-0767. to 8.00 dBresolin et aUl2007h 
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The stellar component of IC 16 13 is known to exten d lit- 
tle farther than 16'. 5 from its centre dBernard et al.L |2007|) . The 
optically brightest part is well separated into 2 sides by a cen- 
tral strip of Hi, usually referred to as 'the bar'. The bar is the 
NE ridge of a hole of the Hi distribut i on wh ere the SW side 
of the galaxy is enclosed (iSilich et"an, l2006l) . Star formation 
is ongoing at the two sides of the bar and at the SW part of 
the Hi cavity, but it is more intense on the NE lobe of the 
galaxy wh ere spectacular giant Hn shells blown by OB as- 
sociations jMeaburn et al.L Il988t IWdez-Gutierrez et all 1200 lb 
are seen. The gal axy has one of the few know n WO of the 
Local Group (W0 3 



Table 1. Log of INT-WFC observations of IC 1613 



Kingsburgh & BarlowljlQQS V discovered by 



iDodorico & Rosal ( 1982 )) and a supernova remnant discovered 
bv lDodorico et al.l ( ll980l) . 

ISandagd (1197 lb presented the first catalogue of targets 
in IC 1613, based on Baade's work on photographic plates. 
Subsequent photometric studies with ever improved detectors 
analysed also the galaxy's population from color-magnitude 
diagrams (Hodge, 1978, 1980; Freedman, 1988; Hodge et al., 
1991; Tikhonov & Galazut dinoval I2002t IBernard et al.L 120071: 
Battinelli et al.. 2007, among others). The DM is well known 



from numerous studies of Cepheids and RR Lyrae (see for in- 
stance 'UdaisH]etal](|200ll))- The youngest stars (age ~ 10 Myr) 
are patchil y arranged on the uniform underlying distribution of 
red stars (lFreedmani[l988l:lBattinem et all 120071) . In the bubble 
region, the young populatio n can be split into 2 a ge groups of 
5-20 Myr and 200-300 Myr (iGeorgiev et aUll999l) . The red su- 
pergiant stars in the region are massive (2O-25M0) and belong to 
the younger age group (Borissoyaet al., 2000). 

Spectral types have been determined for the brightest stars 
of the galaxy and some blue stars in the active bubble region . 
The earliest known star is an 03-04V (Bresolin et al., 1 20071) . 
Inside the H ii super-bubble s , one O f star and three O giants were 
found bv iLozinskava et aTl ( |2002|) . but a larger population of O 
stars is expected: to reproduce the observed far infrared emis- 
sion by dust, it must absorb the equivalent UV emission of about 
9 06-type stars (Meaburn etjiL, 19_88). Recent quantitative spec- 
troscopic anal yses produced the stell ar parameters of a sample of 
B supergiants dBresolin et al.L l2007l) . The upper main sequence 
seems to be well populated but the census of early-type massive 
stars is not yet complete. 

The ultimate goal of our project is to provide a compre- 
hensive characterization of the massive stellar population of 
IC 1613. In preparation, this paper presents a deep photometric 
catalogue (down to V~25) that extends well beyond the visually 
brightest centre of the galaxy, covering a total field of view of 
34' X 34'. The catalogue is provided with a very accurate as- 
trometric calibration, suitable for follow-up multi-object spec- 
troscopy with fiber-based instruments. To understand the pop- 
ulation of blue massive stars, we have built a new list of OB 
associations u sing our automatic friends-of-friends code (after 
lBattineliilll991) . The improvement with respect to previous simi- 
lar works is the input photometric catalogue, that provides a bet- 
ter coverage of the galaxy, reaches a fainter V-magnitude than 
previous ground-based works and has an improved astrometric 
solution. 

The observational dataset and the reduction process are de- 
scribed in Section |2] The astrometric solution is provided in 
Section [3] Our code for automatic search of OB associations is 
described in Section]?] and its application to our IC 1613 cata- 
logue of candidate OB stars in Section|5] The resulting list of OB 
associations is presented in Section|6] In Section|7]we discussed 
how a different input catalogue of stars may alter the results. A 
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comparison with previous catalogues is provided in Section [8] 
Finally, we present our conclusions in Section|9l 



2. OBSERVATIONS AND REDUCTION PROCESS 

The observations were carried out on the 2.5m Isaac Newton 
Telescope (INT) using the Wide Field Camera (WFC). The WFC 
consists of 4 chips (each one is a thinned EEV CCD, with 
2048 X 4096 pixels) that together cover a field of 34' x 34' and 
has a pixel scale of 0.33"/pixel. The log of the observations is 
provided in Table[T] They were structured in 3 sets that included 
300s exposures taken with filters B VRI, and 600s exposures with 
filter U. Each set diff'ers from the other on the pointing coordi- 
nates: second and third are offset by ± 30" in declination so no 
star was lost because of the inter-chip gaps or CCD deffects. An 
additional short exposure of 60s was taken for all the filters in 
the first set, to properly acquire bright sources. 

The S/N achieved for an apparent magnitude of V=23 is: 21 
/ 52 / 37 / 27 / 14 for UBVRI respectively. As we show in the 
next Section, the maximum V-magnitude error for a V=23 star 
is 0.1. The average full- width half-maximum (FWHM) of point 
sources is 3.5 / 3.2 / 2.8 / 2.6/ 2.6 pixels for UBVRI i.e., 0.86"- 
1.16". At the distance of IC 1613, 1" is equivalent to 3.5 pc. 

Complementary VLT-VIMOS images were used for inspec- 
tion purposes. The observations are the pre-imaging part of pro- 
gramme 078.D-0767, Principal Investigator A. Herrero. 



2.1. Photometric Reduction 

The photometric reduction was kindly performed by R Stetson. 
The data were first corrected for non-linearity, following the pre 
scriptions given in the web pagf] of the WFC. Standar d IRA 
routines were then used for bias and flat field correction (IValde; 
19971) . No defringing mask was applied to the I band images 
since no reliable map could be obtained from this dataset. 

The source extraction was performed on the four chips in- 
dependently. The DA OPHOT/ALLFRAME softwai-e package 
(IStetsonLlT987lll994 was used according to the following pro- 
tocol. A first search for stellar sources and aperture photometry 
was performed on each individual exposure. To build the list of 
stars suitable to model the point-spread function (PSF), bright 
and isolated stars with small photometric errors and good shape 
parameter ( sharpness, hereafter #-parm. ) were selected from 
the catalogue of extracted sources. An average of 100 to 200 
stars were picked per chip, spread throughout its surface to prop- 
erly model the PSF spatial variations. 

Allframe then requires two input files: a list of stars and the 
geometric transformation needed to match different frames. The 



http://www.ast.cam.ac.uk/~wfcsur/ 
- IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Fig. 1. Photometric accuracy and goodness of DAOPHOT fit-to-PSF parameters of targets in our IC 1613 catalogue (black dots) as 
a function of magnitude. Sources registered at different chips of the WFC are plotted in different columns. The catalogue's faint 
limit is V=25. V-magnitude errors for V < 23 stars are smaller than 0.1. Grey dots are the input set of targets for the search of OB 
associations. The horizontal lines in the UBV-errors, #- and;t'-parm. plots represent the quality requirements explained in Sect. lS.l.TI 



latter was estimated using daomatch/daomaster. A new input 
list of stars was extracted from the median-stacked image (built 
from all exposures and all filters). AUframe was run once again. 
Adopting the more accurate stellar centroid and magnitude de- 
terminations by allframe enables (i) a better model of the PSF, 



(ii) a refinement of the geometrical transformation, and (iii) a 
better median image. 

As a by-product, allframe creates individual residual images 
after source extraction. Visual inspection of the median stack of 
these, revealed that many stars were missed in the first search. 
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especially in the central region of the galaxy where crowding is 
quite severe. A second search of stars was performed and the 
new targets were added to the input list of stars of allframe, 
which was run a second time with updated PSF and geometrical 
transformation. This loop (allframe - inspection of the median of 
the subtracted images - update of input list) was iterated a third 
and a fourth time for chip 4, that contains the bulk of IC 1613. 

Finally, the photometric calibration was performed using lo- 
cal standards from P. Stetson's catalogued 

The resulting catalogue has almost 85,000 targets. It will 
be publicly available on-line, with an accurate astrometric cal- 
ibration for the target coordinates (see Sect.|3]l. The photometric 
depth and typical accuracy of the catalogue can be seen in Fig.II] 

The expected foreground contamination in the central 15' x 
15' of IC 1613 is -40 stars with V € [16,24] and B-V<0.8 
dRatnatunga & BahcaUl (Il985h from the Milky Way model of 
iBahcall & Soneiral (Il980h ~ This region, covered mainly by 
Chip 4 and marginally with Chips 1 and 3, contains approxi- 
mately 23,500 objects of our catalogue with the same magni- 
tude and color-range. Foreground contamination is therefore ex- 
pected to be negligible in the centre of the galaxy. Chip 2 cov- 
ers the outskirts. It is the least populated and contains only 235 
entries in the listed magnitude and color range. Even though 
a large populat i on of RGB stars is expected in this location 
jBattinelli et al.L |2007|) . Chip 2 is likely dominated by fore- 
ground stars. 



3. ASTROMETRIC CALIBRATION FOR WFC IMAGES 

The importance of an accurate astrometric solution for any kind 
of source catalogue is often underestimated. Firstly, it is es- 
sential for cross-identification with other existing and future 
databases, allowing the combined use of complementary data. 
Secondly, it enables follow-up observations with instruments de- 
manding high position-accuracy, such as multi-object spectro- 
graphs. Additionally, instruments with a large field of view (such 
as the INT- WFC) suffer from large geometric distortions. Our as- 
trometric calibration ensures that the instrumental defects on the 
target positions are totally removed. 

Typically, the transformation from instrumental coordinates 
into celestial coordinates is accomplished by describing the in- 
strument focal plane with a polynomial model. An existing refer- 
ence catalogue is used to determine the parameters of the model, 
i.e. the coefficients of the polynomial. A significant overlap of 
stars between the observations and the reference catalogue is 
needed in each exposure, to ensure that the model is well con- 
strained. The solution is then applied to all (x,y) detections, 
yielding celestial coordinates for all observed objects and not 
only for the reference stars. 

The routine we used for astrometric calibration is fully 
automatic and does not require a previous manual cross- 
identification of the observed stars in the reference catalogue. 
The programme first makes a preliminar cross-identification 
matching geometrical patterns of both input lists of stars (the 
measures and the reference one) in the frame's coordinate sys- 
tem. A linear transformation is calculated and applied to all stars. 
A new cross-identification is made, this time only searching for 
positional matches with a tolerance of 0'.'5. This way, a more 
complete list of common stars is achieved. They are used as refe- 
rence stars to calculate the final solution, a full cubic polynomial 



transformat ion. A more detailed descrip tion of the procedure is 
provided in lCabrera-Lavers et al.l (l2006h . 

In this work, t he reference posit ions were taken from the 
2MASS Catalogue (ICufi-i et al.Ll2003h whose reference system is 
akeady on the International Celestial Reference System (ICRS). 
The 2MASS Catalogue does provide a dense reference list of 
targets with precise positions. Proper motion correction was not 
necessary since our observations and 2MASS' are contempo- 
rary. Since 2MASS is an infrared survey, we cross-correlated 
our catalogue's R magnitude with the J magnitude of 2MASS. 
The starting point was the instrumental (x,y) positions derived 
from the PSF fitting described in Section lTTI 



3.1. Mosaic System and Distortion Removal 

Detectors with field curvature introduce geometric distortion 
into the images. The image will experience magnification to a 
factor that depends on the distance to the optical axis. The result- 
ing radial displacement is usually proportional to r^, although 
perfect azimuthal symmetry is not always seen. This eff'ect must 
be corrected prior to any celestial coordinate calculation. 

It is often found in the literature that the radial distortion of 
the INT- WFC is removed by following the instructions given in 
the instrument web page. It suggests performing a radial cor- 
rection with origin in the optical center of the mosaic (i.e., the 
composition of the four chips) to calculate a polynomial solution 
afterwards. 

We first transformed the instrumental (x,y) coordinates of all 
targets on all CCD into global mosaic system (Xn,, ym) coordi- 
nates with origin in the optical center, using the linear equations: 



x„ — a*x + b*y + c 



(1) 




■am 



' http://www3.cadc-ccda.hia-iha.nrc- 
cnrc.gc.ca/community/STETSON/standards/ 



x(pix) 

Fig. 2. Stacked position residuals of the complete INT- WFC mo- 
saic. The arrows represent the difference between the positions 
listed by 2MASS and our catalogue, for the reference stars. The 
scale vector of l"size is plotted on the lower right corner A ra- 
dial correction plus a global cubic polynomial solution has been 
applied to our catalogue. The positional differences vary as a 
function of the mosaic coordinates x and y. The fixed pattern 
shows that instrumental distortions have not been completely re- 
moved. 
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Fig. 3. Stacked position residuals as a function of the chip coordinates (x,y) after radial correction. The (J! 1 size vector is plotted on 
the lower right corner of the panels. Only a linear characterization of the detector is calculated and applied. The fixed pattern shows 
that an additional correction is still needed, and that the correction is different for each CCD. Similar plots for Chip 1 and Chip 2 
present different distortions. 



where the abode f coefficients for each CCD were taken from 
the technical reports on the WFC. 

The radial correction to the coordinates (Xm, ym) is defined 
as follows: 



Ax^R x,„ (xl + yl) 



Ay^Ry,„{xl+yl) 



(2) 



where R=-5.30 x lO '^pixel"^ is the distortion coefficient. 

After the radial correction was applied, we used a general 
cubic polynomial solution to describe the detector, as explained 
in the previous Section. We chose a third-degree polynomial so 
additional distortions could be accounted for. 

The residuals of the reference stars are the difference bet- 
ween the transformed instrumental positions and the reference 
catalogue coordinates. In Figure |2] we have plotted the vector 
field of residuals as a function of position in the focal plane. A 
fixed pattern geometric distortion is still present in the INT- WFC 
system. After several failed tests, we concluded that the com- 
bination of a global mosaic radial correction and a subsequent 
general polynomial solution cannot fully correct the distortion 
caused by this telescope-camera combination. 

The reason why our characterization failed is that each chip 
demands a different correction. Figure[3]shows that after correct- 
ing for radial distortion systematic errors remain in the relative 
positions (as expected) but each CCD presents a very different 
pattern of distortion. Hence, for a more accurate calibration we 
deviated from the generally adopted strategy: after correcting for 
radial distortion, we estimated an individual fourth-degree poly- 
nomial solution for each chip as opposed to calculating a global 
solution for the whole mosaic system. 

3.2. Results and Error Estimate 

We estimate the uncertainty of our calibration by comparing our 
final catalogue with a high precision one: 2MASS. We plot the 
difference between our positions and 2MASS' in FigurelH as a 
function of the R-magnitude of the targets. The difference be- 
tween both catalogues is larger in fainter objects but the average 



is zero, as expected. The standard deviation of the difference is a 
more meaningful parameter If the complete catalogue is consid- 
ered, the standard deviation is (cTacose, erg) = (CK'IQ, 0'.'16). Note 
that the standard deviation is larger ((X'S) when a global solu- 
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Fig. 4. Difference between our catalogue and 2MASS positions 
for the reference stars, as a function of R-magnitude. Upper 
panel: difference in right ascension. Lower panel: declination 
differences. The dark line represents the standard deviation. 
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tion for the whole mosaic is calculated. The number of reference 
stars in each CCD is 140, 123, 138 and 206 stars, respectively. 

Let us also remark the lack of any magnitude trend of the 
differences: the standard deviation increases for faint objects, 
but the average is always ~0. This fact evidences that the po- 
sitions obtained are free of any obvious "magnitude equation", 
the undesired correlation between the target position relative to 
the image center and the magnitude of the target. We also inves- 
tigated the dependence of position difference versus color, and 
no systematic difference was found. 



4. AUTOMATIC SEARCH OF OB ASSOCIATIONS 

Catalogues of extragalactic OB associations meant as clustering 
of bright blue stars over a haze of field stars can already be found 
in the first extragalactic photographic jobs. However, marking 
OB associations 'by eye' is a subjective procedure; the results, 
namely the sizes of the found associations, depended on a num- 
ber of factors: distance to the galaxy, p late sca l e, exp osure time 
and selection criteria (see discussion of lHoded (Il986h ). 

The advent of CCDs combined to target extraction rou- 
tines (such as DAOPHOT) has given rise to the automatic 
search of OB associations. The algorithms usually follow two 
main philosophies: (i) search for statistically si gnificant peaks 
(over 2 or 3 cr) in m aps of stellar density dKontizas et al.L 



Table 2. Photometric-quality requirements for stars to enter the 
algorithm. 



1199 6; Gouliermis et al., 2003) or (ii) the Path Linkage Criterion 
Method (Battinelll, 1991) also known as friends-of-friends algo- 
rithm . The use of pattern recognition dFeitzinger & Braunsfurthl 
11984 is less common. 

We have developed an algorithm based on the friends-of- 
friends method: two stars belong to the same association if the 
distance of one star to the next is less than a given amount. The 
input entries to the code are a list of points with spatial coordi- 
nates and two free parameters, the search distance (Ds) and the 
minimum number of members (Nmi,,) for a group to be consid- 
ered an association. 

Given an OB star, the programme looks for other OB tar- 
gets within the search distance and registers them as members of 
the same association. The search is repeated for all the new OB 
members until no star is found within Ds from any of the periph- 
eral stars. In the process, non-OB stars are also included in the 
association if they are less than Ds away from an OB member 
Groups with less than Nj^in OB members are discarded. 

The one case the algorithm would not detect is a chain of 
stars each closer to the next less than Ds forming a ring of ra- 
dius larger than Ds. The code would walk the ring and miss the 
central targets. At the moment, this has to be checked visually. 

There is no physical constraint for the parameters Ds and 
Nmin- In fact, friends-of-friends users do not agree on how to de- 
cide their values. Nmin is more problematic, since the expected 
population of an as sociation depends on the mass of its parent 
cloud and the IMF. iBorissova et alj (2004) decided Nmin upon 
evaluating the chances that an association is actually an acciden- 
tal group of stars. Other works consider that an OB association 
is valid if its stellar surface density (E) is statistically ab ove the 
mean, and choose Ds and N^in on these grounds. Ivanovl ([19961) 
uses Ds that maximizes the mean normalized fluctuation of E. 
IWilsonI ( Il991i) sets Ds = 1 / VTtt x S) and N^in by compari- 
son of the resu l ting a ssociations to a random area of the galaxy. 
iGeorgiev et al.l d 19991) find Ds and Nmin from the Fourier trans- 
form of the surface stellar density. 

In this work we infer Nmin from Milky Way associations, 
by counting the members observable if they were located at the 
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dist ance of the studie d gal axy (see Sect. | 5.2|). For Ds , we fol- 
low iBattinellil (Il99lh and iBattinelh & DemersI (1 19921) and set 
this parameter to the abscissa of the maximum of the function 
fy,„,„{Ds), the number of associations that the algorithm finds 
with at least Nmin stars when run with search distance Ds. At 
larger Ds associations tend to merge into a single larger system. 
On the other hand, if Ds is small then the associations separate 
into their individual components (isolated stars or binaries). We 
adopt this criterion because it is reproducible, but we stress that it 
does not guarantee that the groups found in this way correspond 
to physical associations. Note that the visual identifications of 
OB associations already encountered this problem, and deciding 
whether a star clustering is only one association, or it can be split 
into two or more smaller ones was quite arbitrary (Hodge, 1978!). 
iHodgel (119861) already pointed out that additional information is 
required to make sure that we are dealing with physical groups. 
In fact, the same problem is found when defining OB associa- 
tions in the Milky Way. 

Bastian et al. (2007, |2009|) have recently shown that there 
is a correlation of the typical size of the associations found by 
friends-of-friends algorithms and the Ds value used. These re- 
sults stress the fact that the OB associations identified this way 
are not necessarily bound physical entities. Nonetheless, enclos- 
ing OB stars in groups is a useful and necessary first step to 
study the youngest galactic population. Further work studying 
the color-magnitude diagrams of the associations found in this 
work (Garcia et al., in prep.) will allow us to progress towards 
the identification of the real associations. 



5. OB ASSOCIATIONS IN IC1613 

The catalogue of OB associations was made in 3 steps. We firstly 
built an input catalogue of candidate OB stars (Sect. lSTI i. Values 
for the two free parameters of the code, Ds and Nmin, were de- 
cided afterwards (Sect. 15.21 ). Finally, the code was run to pro- 
duce the list of associations, which we inspected individually 
(Sect.O- 



5. 1 . Input Catalogue of Stars 

The starting point is the deep photometric catalogue described in 
Sect. 121 We kept only stars with exceptional photometric quality, 
setting hard constraints on the errors of U, B and V magnitudes. 
The reason is that we built the list of likely OB stars in IC 1613 
based only on their reddening-free Q color (see below), a linear 
combination of U, B and V. If allowed errors in these magnitudes 
were as large as 0. 1, the resulting Q error would be 0.2 and could 
lead to target mis-classification. We required errors smaller than 
0.05 magnitudes. Note that by setting hard constraints on UBV, 
we do not include the red population with faint U magnitudes 
(and large U-errors). 

We set looser constraints on the DAOPHOT parameters that 
describe the goodness of the PSF extraction of each target. We 
allowed for a #-parm. interval (roundness of the PSF) of [-2, 2]. 
The X parameter (quality of the PSF fit, hereafter ;i'-parm.) in- 
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Fig. 5. Locus of O and B stars in the U-B vs Q diagram. 
For each Q-interval delimited by a horizontal band, the cen- 
tral box corresponds to unreddened or moderately reddened ob- 
jects. The symbols re present stars of different spectral type (from 
iBresolin et al.l (l2007h ) as detailed in the legend. Most O stars and 
B supergiants are found in the central unreddened box of the two 
bottom strips. 



creases with magnitude and tight constraints on this parame- 
ter may leave out the brightest stars. In fact, two trends of the 
X-paim. are seen in Fig. [1] corresponding to the long and the 
short exposures. DSS, GALEX and/or VIMOS images reveal 
that most of the faint targets with;t' > 5 are background galaxies 
or blends. A very small fraction corresponds to stars that satu- 
rate. We will allow for;^' ^ 4 for faint stars, and;^' < 30 for bright 
targets. 

The photometric qualification criteria are summarised in 
Table |2] Fig. [T] shows the typical photometric errors of the com- 
plete catalogue (black), with the targets that comply with the 
photometric requirements marked in grey. From a total of 84,771 
targets in the initial catalogue, we keep 3,960. 

5.1.1. Candidate OB Stars 

Candidate blue massive stars were chosen from their Q color, the 
reddening-free parameter Q - (U - B) - 0.72 x (B - V). Q in- 
creases monotonically towards later spectral types in the interval 
Q 6 [-10, -0.4], corresponding to O3-A0 supergiants and 03- 
B5 dwarfs (see for instance the calibration of Fitzgerald ( 1970)). 
We defined the locus of O and B stars in the U-B vs Q dia- 
gram (Fig. [STi using the co l ors of stars with spectral types de- 
termined bv lBresolin et alj ( l2007l) . While Q is an indicator of 
spectral type, U-B holds the information of whether the star is 
reddened or not. Most O and B supergiants are found in unred- 
dened boxes with Q < -0.4. IC 1613's known WO is found in 
the diagram as an outlier in the lower-left corner 

Q was the only selection criteria to find OB stars. We allowed 
for wide ranges of U-B and B-V to account for reddening since 
we suspect that IC 1613's internal extinction is not negligible, 
particularly for OB associations, and that the dust content may 
be significant in the bar (see full discussion in Sect. 15.1.21 ). 



Table 3. Photometric constraints for OB stars 



Parameter 


MIN 


MAX 


V 


16 


24 


Q-parm 


-2.0 


-0.4 


U-B 


-2.0 


2.0 


B-V 


-0.5 


2.0 



Table 4. Photometric constraints for non-OB stars 



Parameter 


MIN 


MAX 


V 


16 


24 


Q-parm 


-0.4 


1.0 


U-B 


-2.0 


2.0 


B-V 


-0.5 


2.0 



The allowed V interval is V e [16,24]. This range ensures 
that all the interesting targets for our spectrosc opic studies ar e 
included: 03 I/BO V types, with Mv=-6.5/-3.8 ( lMassevi[T998h . 
At the adopted distance modulus DM=24.27, the observed mag- 
nitudes will be V= 17.8/20.5. We allow for extra 1.5 magnitudes 
in the bright interval limit, not to leave out exceptional cases 
such as hypergiants. The limit for faint stars is also extended to 
allow for reddening and include dwarfs of later B-types. 

Table [3] summarizes the photometric properties of candidate 
OB stars. The photometric constraints set on non-OB stars are 
listed in Table |4] This is not an attempt to determine the ex- 
act spectral sub-types of OB stars from Q. The degeneracy of 
blue massive star colors at optical wavelengths may be broken 
in the future with improved phot ometric capabilities c ombined 
with new statistical techniques ( Mafz-Apellanizl |2004|) . but the 
required accuracy (one thousandth of a magnitude) is not feasi- 
ble to date. 

We note here that the use of the Q parameter has two caveats, 
(i) It is reddening-independent as long as the reddening law is 
standard, (ii) Beyond very early spectral types the univocal be- 
haviour of Q cannot be ensured. At Q — 0.4, it varies rapidly with 
spectral type and luminosity class. The inherent scatter of the 
calibration of Q with spectral type may take mid/late B dwarfs 
and late-B/early-A supergiants to values larger than -0.4. This is 
the case for the three "late-B" stars in the Q>-0.4 box of Fig.|5] 
Late stars around G-type are close to a local minimum of the 
Q vs spectral type relationship, at Q — 0.37. The scatter of the 
calibration or specific stellar properties (UV excess from stellar 
activity, for ins tance) may bring the s tars to the boxes of the early 
types. Finallv. iBresoUn et all J2007h 's A8 (A2 la) is also in the 
Q<-0.4 locus of B stars. The authors found H„ and H/j P Cygni 
profiles, and their published spectra shows that the Balmer lines 
in the blue and near-UV are unusually weak. Clearly, this star is 
peculiar and hence not expected to follow the general trend. 

We find that the candidate OB stars concentrate in Chip 4, 
mostly on the NE lobe of the galaxy where the giant bubbles 
are seen, but also in the opposite direction, in the SW side of 
the giant Hi cavity where part of the galaxy is enclosed. This 
is in agreement with previous results (see Sect. [TJ- Despite our 
large allowance for V, U-B and B-V to include the most reddened 
targets, we find very few blue stars in the bar 

The 20 brightest blue stars from our catalogue of OB candi- 
dates are listed in Table |5] Spectral types are known for 12 of 
them. This subset includes an A dwarf and a G2 star (the two 
"A and later" stars w ith Q<-0.4 in Fig.[5]l, which are foreground 
objects according to IBresolin et al.l ( l2007h . It is remarkable that 
only one red star has entered Table |5] supporting our color se- 
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Table 5. Brightest candidate OB-stars: magnitude and colors. Column 'ASSOC indicates the associatio n the stars belong to a s found 
in this work (-1 means they are isolated). We also provide identification numers and spectral types from lBresoUn et all (|2007|) . when 
available. 



ID 


RA[deg] 


DEC[deg] 


V 


B-V 


U-B 


Q 


ASSOC 


Bresolin et al. (2007) 




(J2000.0) 


(J2000.0) 












ID 


SpType 


39916 


16.194780 


2.166142 


16.43 


0.54 


-0.05 


-0.44 


-1 


C12 


G2 


59802 


16.246347 


2.154164 


16.44 


0.12 


-0.58 


-0.67 


127 


A8 


A2 la 


83197 


16.369705 


2.083662 


16.67 


0.19 


-0.51 


-0.64 


-1 


_ 


_ 


52033 


16.224657 


2.119503 


17.24 


0.63 


-0.02 


-0.47 


-1 


_ 


_ 


62390 


16.253813 


2.178189 


17.42 


-0.14 


-0.98 


-0.88 


137 


AlO 


Bl la 


69336 


16.276543 


2.158706 


17.69 


-0.15 


-0.98 


-0.87 


185 


B3 


BO la 


15779 


16.124454 


2.051565 


17.76 


0.64 


-0.10 


-0.56 


-1 




_ 


42294 


16.200409 


2.079537 


17.95 


0.65 


0.01 


-0.46 


60 




_ 


65875 


16.264278 


2.162671 


18.06 


0.01 


-0.55 


-0.55 


161 


\ 1 


B9Ia 


1705 


16.023046 


2.209241 


18.13 


0.44 


-0.16 


-0.48 


-1 






60449 


16.248133 


2.154477 


18.23 


-0.13 


-0.97 


-0.88 


127 


B4 


B1.5 la 


56665 


16.237545 


2.084683 


18.24 


0.37 


-0.41 


-0.68 


114 






23263 


16.151031 


2.190748 


18.24 


0.39 


-0.24 


-0.52 


-1 


C5 


A V 


58040 


16.241360 


2.096872 


18.39 


-0.09 


-0.88 


-0.82 


120 


A7 


B2 lab 


27381 


16.163216 


2.146994 


18.42 


-0.07 


-0.88 


-0.83 


-1 






68599 


16.273815 


2.096740 


18.54 


-0.09 


-0.85 


-0.78 


-1 


B13 


B5 lab 


65133 


16.261927 


2.083119 


18.56 


-0.05 


-0.70 


-0.66 


163 






52841 


16.226801 


2.136780 


18.62 


-0.07 


-0.72 


-0.67 


105 


A5 


B3Ib 


34964 


16.182615 


2.112548 


18.62 


-0.13 


-1.04 


-0.94 


43 


BU 


09.5I 


18570 


16.135248 


2.158957 


18.68 


-0.13 


-0.94 


-0.85 


21 


C6 


Bl la 



lection criteria. Despite its red intrinsic colors the G2 star is in 
the limit of the accepted Q interval, justifying the need of hard 
constraints on the photometric errors of stars. The remaining 10 
targets (3 of them brighter than V<17.8) are O, B and A super- 
giants, with high chances of been IC 1613 members. Table |5] 
supports our selection method and also indicates that the fore- 
ground contamination is small. 

5.1.2. Variable Extinction and Extinction Map 

The foreground extinction towards IC 1613 is remarkably low 
(E(B-V)~0.2, see Sect.[T]). In the past it was argued that the in- 
ternal reddening was also negligible since background galax- 
ies can be seen throug h IC 1613 (Hodge, 1978; Fre edmanl 
ll988HKrienke & Hodgd 12000). although Hodge (1978.) already 
marked the position of ~10 possible dust clouds in the galaxy. 

More recent papers yield higher val ues of ex tinction in 
the lines of sight of individual objects. [Humphreys (1980) 
derived E(B-V)=0.05-0.10 towards early type members, and 
~0.2 towards M- supergiants. The ceph eid studies of the OGLE 
dWillick & Batrai r2001) and Araucaria dPietrzvnski et alll2006l) 
projects yield average extinction values towards their targets of 
<E(B-V)>= 0.085 ±0.052 and 0.090±0.019 , with extinction law 
similar to the MW. iTautvaisiene et al.l(l2007l) 's study of three M- 
supergiants yields three different val ues E(B-V)-0 .04, 0.06 and 
0.03. In the bubble region, G eorgie v et al.l (1 19991) report E(B- 
V)=0.06 and Lozinskava et al.l(l2002[) finds E(B-V) in the range 
of 0.18-0.4 towards 9 OBA stars , considerably la rger The ex- 
tinction towards the WO is 0.1 (iLee et all l2003h . but redden- 
ing is suspected to be high and anomalous towards Wolf-Rayets. 
Detections of the molecular cloud, neutral hydrogen and dust by 
IRAS and Sp i tzer in the bubble region (iMeaburnetaLL 119881: 
I Jackson et al.L 12006) are further evidence of extinction in the 
area. 

We therefore suspect that internal reddening is highly vari- 
able and far from negligible. To prove it, we have estimated the 
color excess in the lines of sight of our candidate OB stars in 
IC 1613. Intrinsic colors were calculated from Q, thanks to its 



univocal relation with spectral type (hence (B - V)o) that holds 
for the early-type s tars. This relat i on wa s parameterized for low- 
metalUcity stars bv lMassev et all (l2000h : 



(B-y)o = -0.005 H- 0.317 X e 



(3) 



using Kurucz' ATLAS9 2=0.082© models. (No such relation- 
ship exists for metallicity as low as IC 1613's, as far as we 
know). This is an approximation, since equation [3] does not 
hold if the reddening law deviates from the standard behaviour. 
Another limitation of equation [3] is that it is independent of lu- 
minosity class. Although convenient for fast calculations, it will 
yield (B - V)() slightly bluer than the actual intrinsic colors of 
supergiants and slightly redder for dwarfs. The effect will have 




Fig. 6. Extinction distribution towards the blue massive star can- 
didates of IC 1613. Foreground reddening has not been cor- 
rected. Most of the stars have E(B-V) in the range 0.05-0.25 and 
the peak of the distribution is located at 0.10-0.15, significantly 
higher than the published foreground extinction. 
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Fig. 7. Extinction maps. Only the portion of IC 1613 containing OB associations is shown (same as in Fig. |9]l. Different panels 
consider different E(B-V) ranges. Clock-wise starting from top-left: (1) complete range O.-l.O, (2) the traditionally considered 
range for this galaxy O.-O.l, (3) intermediate E(B-V) range 0.1-0.5, and (4) large extinction 0.5-1.0. 



a maximum for B-type stars. For instance, th e comparison of 
the E(B-V) provided bv lBresolin et all ( |2007|) for seven early- 
B supergiants (0.06 on average) with the values derived from 
equation |3] (0.1 5 on average), hints that the latter overestimates 
extinction for these stars. However, we expect that statistically 
these effects will compensate in the large stellar population of 
the galaxy. 

The distribution of the extinction towards the blue objects 
of IC 1613 is shown in Fig. |6] It peaks in the E(B-V) interval 
0.10-0.15. Most of the stars have color excess ranging 0.05 to 
0.25, higher than the foreground extinction provided in the lit- 
erature. For a few objects we obtain E(B-V)>1. Most of them 
(more than 90%) are background galaxies for which equation [3] 
does not apply, and have been discarded. The negative values of 
E(B-V) in Fig. |6] are due to a combination of factors, mainly 
the intrinsic scatter of the calibration of the Q parameter with 



(B - V)(), but also local changes in the slope of the reddening 
law or photometric errors. 

The spatial distribution of E(B-V) is shown in Fig. [T] The 
different panels are the extinction maps for 4 different extinction 
intervals: the complete range (O.-l.O), the traditionally consid- 
ered range for this galaxy (O.-O.l), intermediate interval (0.1- 
0.5), and large extinction (0.5-1.0). Starting from the pixel map 
of the INT-WFC Chip 4 image, stars were distributed into boxes 
of 25 X 25 pixels. Each box was assigned the average extinction 
of the enclosed stars with E(B-V) in the studied range. Empty 
boxes (i.e., bins with no stars with E(B-V) in the considered 
interval) where assigned the minimum of the interval. In the O.- 
l.O and O.-O.l panels, empty boxes are assigned the minimum 
positive value found in the sample E(B-V)=0.0002. 

Boxes with E(B-V) in the range O.-O.l trace well the star 
forming regions of the galaxy, which is also the preferred lo- 
cation of OB star candidates. A concentration of E(B-V)~0.1 
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Fig. 8. Number of OB associations found as a function of Ds and 
Nmin, fi\[„,i,XDs)- Note the irregular scale on the x-axis. Different 
symbols represent the results for different values of N^in as in- 
dicated in the legend. There is a local maximum around Ds=6". 
The slight increase at Ds=50" is caused by stars in the galac- 
tic outskirts, far away between them, that cluster because of the 
large Ds admitted value. 



points is seen at the bubbles and at the SE edge of the galaxy. 
Points with larger extinction values (0. 1 <E(B- V)<0.5) also trace 
the star forming regions, with a slight concentration towards the 
bar area and the SW part of the H i cavity. Larger E(B-V) values 
are distributed almost randomly throughout the galaxy. These 
results support that extinction is not uniform, with rapid spatial 
variations. 

5.2. Free Parameters of the Code 

We followed the procedure described in Sect. |4] to determine 
Ds and Nn^n^ 

We used l Humphreysl ( ll978l) 's catalogue of Galactic OB as- 
sociations to estimate N^in. We removed all luminosity class- 
V stars with spectral type later than B5. ( According to the 
V-magnitude limits of Table [3] the catalogue of OB candidates 
would include spectral types down to B8 V, but the Q constraints 
set the cut at B5 V). 25 out of 78 associations (i.e., ~3Q%) are 
left with 3 or less OB members. Since projection effects may be 
at work, we consider that 2 members are too few for an extra- 
galactic OB association. We therefore set the minimum number 
of stars to 3, at the risk of neglecting the less massive associa- 
tions. We expect that we will be able to spot out false OB associ- 
ations through inspection of the images, and the color-color and 
color-magnitude diagrams. 

We ran the code for different search radii in order to as- 
sess /Af„„„(-Ds) (see Fig. [S), and set Ds to the abscissa of the 
maximum. Besides Nniin=3, we tested three additional values: 
Nmin=4, 5 and 6. The bump of fy^^.^X^s) flattens as N^in in- 
creases, but, in the four cases there is a local maximum at Ds=6- 
7". This local maximum is preserved up to Nn,in=10 (not shown 
in Fig. [8]). 

We therefore set Ds - 6" and N^in = 3. 

Ds=6" is similar to the values derived by Borissova et al. 
(IC 1613, 6.6") and Pietrzynski et al. (NGC 300, 7.7"), and 
much larger than those adopted by Bresolin et al. (1996, 1998, 
0.83 to 3.9" in different spiral galaxies at ~10 Mpc). Since 
Borissova et al., Pietrzynski et al. and this work have used 2m- 
class ground-based telescopes with comparable instrumentation 



whereas Bresolin et al. used HST-WFPC2 data, it is tempting to 
interpret the differences in terms of the superior plate scale of 
the WFPC2 (O.r'/px) breaking the OB associations into smaller 
pieces. However, the physical values of Ds vary between the 
21 and 23 pc obtained by Borissova et al. and ourselves for IC 
1613 and the 35-84 pc obtained for spiral galaxies by Bresolin 
et al. and Pietrzynski et al. Therefore, the smaller angular value 
of Ds obtained by Bresolin et al. is a consequence of the longer 
distances towards their sample galaxies (hence smaller angular 
diameter of associations). Whether the smaller physical values 
found for IC 1613 are due to inherent differences of star forming 
regions of irregular galaxies versus spirals, the distance to the 
galaxies or to poor statist i cs mu st be clarified in future works. In 
the light of lBastian et al.l (l2007[) 's results we also anticipate that 
the different Ds values used will result in different typical sizes 
for the associations, but we leave the discussion on whether or 
not there exists a universal length-scale of OB associations for a 
subsequent paper 

5.3. Final OB Association List 

We ran our automatic search algorithm on the catalogue of can- 
didate OB stars one last time with Ds - 6" and Nmi„ - 3, and 
produced the list of OB associations presented in Sect. |6] All 
stellar groups provided by the code were accepted as valid asso- 
ciations, as long as they have at least 3 valid members (i.e., not 
background galaxies, blends or bright red objects with Q — 0.4). 
Each association was studied thoroughly using images, color- 
color (B-V vs U-B and U-B vs Q) and color-magnitude (B-V vs 
V and Q vs V) diagrams to discard fake detections. The analysis 
of these diagrams will be presented in a second paper. 

We examined high resolution images of IC 1613 from 
the HST archive and VLT-VIMOS using the Aladin software 
jBonnarel et al.L |2000l) to discard non-valid association mem- 
bers, find neighboring associations and/or nearby objects that 
did not enter the input catalogue because of the quality of their 
detection (poorer than the requirements of Table |2]l. 

We also cross-matched our catalogue with published spectral 
classification for IC 1613 targets. The known spectral types of 
association members are provided in Sect. 16.1! 

The code finds 198 groups. 5 of these groups are galaxies. 
After visual inspection, 29 associations were discarded because 
they did not reach the minimum number of valid members and 
47 were marked dubious because we could not assure that they 
had at least Nmin isolated blue stars. The remaining 1 17 groups 
(~60%) are bona-fide OB associations. 



6. CATALOGUE OF OB ASSOCIATIONS 

The list of OB associations built with our code is presented 
in Table IA.13I Their positions are marked on the galaxy in 
Fig. |9] Zooms into different galactic sectors are provided in 
Appendix lA.3l 

The associations with more than 10 OB members are listed in 
Table|6] as illustration of the complete catalogue. The content of 
different columns is: (1) OB association identification number 
(2) and (3) Right Ascension (RA) and Declination (DEC) of the 
centre determined from the average of the positions of the OB 
members, equinox J2000.0 (4) Number of OB members (Nqb) of 
the association. (5) Total number of members registered (Ntot); 
note that our photometric selection criteria is biased against red 
objects because we demand high accuracy in U magnitudes. (6) 
Quadrant allocation in Fig.|9] (7) Notes: 'Gal' - the association 
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Table 6. OB associations with more than 10 OB-candidate mem- 
bers. 



ID 


RA [deg] 


DEC [deg] 


Nob 


Ntot Quad 


Notes 




(J2000.0) 


(J2000.0) 








(1) 


(2) 


(3) 


(4) 


(5) (6) 


(7) 



n 
n 
n 
n 



n 
n 
n 
n 
n 



n 
n 



n 

n 
n 



is actually a galaxy; 'd' - association discarded because it does 
not have enough valid members (i.e. visual inspection reveals 
that some are blends or gas knots); '?' - dubious association 
(associations with a small number of members, with a significant 
fraction under suspicion of being blends or not blue stars); 'n' - 
an entry in Sect. 16.11 

The code also provides us with a list of the OB members of 
each association and also the less massive population. Defining 
OB associations by their members, rather than by hand-drawn 
lines on photographic plates, allows a direct comparison of re- 
sults between different works and with observations. The mem- 
ber listings will be available on-line; an example is provided in 
Table |A:fil 

6.1. Comments on Individual Associations 

Assoc. #10: an A3 II dBresohn et al.l(l2007h 's C3) is located 6".2 
away from one association member 

Assoc. #16: Part of the association shows radial distribution 
around one red and two blue central stars, that were not included 
in the OB catalogue because of their photometric errors. 

Assoc. #21: The brightest blue star is a B 1 la who s e spec trum 
displays strong H„ wind emission jBresolin et al.l (l2007h . id: 
C6). 

Assocs. #27 and #28: their members are horizontally aligned 
in our images (they have very similar RA), perpendicularly to 
a spike of a nearby saturating bright star The brightest stars of 
both associations fall right on the spike, and we discard them 
because their photometry may be a ltered. 

Assoc. #43: The brightest star is lBresohn eTaLl(l2007l) 's Bll, 
an 09.5 I. 

Assoc. #48: Blue stars 7" away in the SE direction could 
connect this ass ociation with #5 L 

Assoc. #53: 'Breso hn et all (l2007 Vs C16, an A5 lb, is in- 
cluded as a red member. 

Assoc. #54: A bright blue nearby (8" away) star is missed 
for the association. 



Assoc. #61: Its brightest member is iBresolin et al.l (l2007l) 's 
B8 (B8 lb). 

Assoc. #69: The brightest star is IBresolin et al.1 (l2007h 's BIO 
(B8 lab). 

Assoc. #75 includes an A5 II dBresohn et al.1 d2007l) 's A2). A 
nearby red star (not listed under the association) could be a red 
supergiant (RSG). 

Ass oc. #83' s bright est blue member has very blue Q. It is a 
Be star dBresohn et al.l ((2007)'s B19); its emission lines may ex- 
plain its bluer color Additional faint blue stars are found nearby 
(at ~ 7"-9"). 

Assoc. #85: An additional blue star is only 7" away. 

Assoc. #86: Its brightest member looks triple in the WFPC2 
images. A nearby bright red star could be an RSG star of this 
association. 

Assoc. #96: iBresohn et all d2007b 's A3 (A7-F0 lb) is only 
7"away. 

Assoc. #97: iBresohn et all (l2007l) 's A4 (B9 II) is only 
9"away. 

Assoc. #98: An additional blue target 9" away is missed for 
the association. 

Assoc. #99: A member with very blue Q has the same posi- 
tion in the U-B vs Q diagram as IC 1613's WO. However, this 
object may be bl ended with th e brightest association member 

Assoc. #103: Breso lin et al.1 (I2007ii) 's A4 star, a B9 II, is in- 
cluded in the association. Because of its late spectral type, this 
object's Q color is slightly beyond our Q constraint for blue stars 
(i.e., Q>-0.4). 

Assoc. #104 contains IBresolin etakl d2007l) 's CI 7 star, an 
early-A II. The target enters the catalogue as a red star (Q >- 
0.4). 

Assoc. #105 contains lBresohn et al.l(l2007l) 's A5 star (B3 lb). 
(This association was discarded because it only has 2 valid mem- 
bers). 

Assoc. #106: The brightest star is IBresolin et al] d2007l) 's 
A4b star (B5 lb). 

Assoc. #114 is a round association located at the outer part 
of the galaxy. The brightest blue star (V~18.5) could be blended 
with the next two brightest stars of the association, although tar- 
gets are well resolved in VLT-VIMOS images. 

Assoc. #116: A BO lab star (.Bresolin el: al.1 (l2007l) 's B14) is 
close to the association. 

Assoc. #117: The position of the red member in the color- 
magnitude diagrams matches that of an RSG. Close to associa- 
tion #114, together they could form a ring. 

As soc. #120: The brightest member is IBresolin et al.l 
(l2007l) 's A7, B2 lab. 

As soc. #122: The brightest blue member is IBresolin et al.l 
(|2007|)'sB16, B1.5 lab. 

As soc. #124's brightest member is an AO III (IBresolin et al.l 
(12007'), id:B15). 

Assoc. #125 : Its brightest member could be contaminated by 
a nearby red star, although its normal colors suggest otherwise. 

Assoc. #126: The brightest star has very blue Q and its col- 
ors deviate from the reddening law, but visual inspection reveals 
that it is a normal isolated st ar. This star i s [AM85]3, with spec- 
tral type OB+em ( Azzopardi et al.Lri988h . The nebular emission 
could explain its anomalously blue Q. 

Assoc. #127 delineates the bright part of one of the bub- 
bles of the NE lobe. #127 may form a single giant association 
with #147. Spectral types are known for three stars in the as- 
s ociation: the brighte st blue stars have types A2 la and B1.5 la 
(iBresolin et al.ld2007l) 's id's are A8 and B4). Skiff (2007) assigns 
AO Ia(e) to the former and additionally lists an RSG, [S71]V32 
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Fig. 9. INT-WFC U-band image of IC 1613 (only the portion of Chip 4 that contains the associations is shown). The positions of 
the associations found in this work are marked in black. The grid divides the image in 9 portions. A zoom into each quadrant is 
provided in Appendix lA.3l Note how the OB associations cluster in sectors 5, 7 and 8. 



with type Ml la (not included in our catalogue because of the 
hard constraints on blue photometric errors). The A-supergiant 
has b een sug gested being a foreground star (Sandage & Katem, 
11976 ). but B resolin et al.l (12007) found its radial velocity in 
agreement with IC 1613, and report P Cygni profiles of Hy 
and H/i (see Sect. 15.1. It . Spectroscopy of higher resolution than 
iBresolin et aLl(l200 7)'s VLT-F0RS2 data is needed to determine 
whether this object is an LBV-candidate. 

As soc. #128: The brightest star is a Be jBresohn et al.l 
(l2007h 's B5), and this may explain its very blue Q color. 

Assoc. #134: A possible additional member, bright and blue, 
is only 7"away. 

Assoc. #136 marks a bubble rim. 



Assoc. #137 is formed by several arc-like structures that pos- 
sibly trace the edges of several bubbles. Four members at the 
West could form a separate association. The b rightest member 
is a B superg iant: Bl la (Bresolin et al.' (2007)'s AlO), B2/3 le 
(ISkiffl, I2OO7I), B2-3 I (Humphrevs, .1980). The second brightest 
star is lBresohn et"ai] (l2007h 's A9 (B5 lab). 

Assoc. #140: A nearby blue star (~7" away) is not included 
in the association. 

Assoc. #144: The brightest member ionizes its surrounding 
nebula. 



As soc. #146: The brightest member is a Be (IBresolin et al.l 
dlOOTh . id: B12). 
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Assoc. #147 is dominated by 2 clumps, each in the inter- 
section of 2 bubbles, plus 2 arc-like structures that trace the 
rim of one bubble each. Some detections are actually neb- 
ular knots (discarded). The association has 5 bright mem- 
bers (V £ [19,20]), four o f them in the largest clump of 
stars. iLozinskava et al.l (|2002|) obtained low resolution spectra 
{AA =8A) of the stars in the main clump, with a 6m telescope 
(note that this area is hardly resolved, even in the VLT-VIMOS 
images). The group includes an Of star, two O giants and two 
OB candidates. The association has an additio nal OB member, 
far fr om the regions of crowding (OB+em?, lAzzopardi et al.l 

dMi). 

Assoc. #148: A nearby blue star, potential member of the 
association, is 7.5"away. 

Assoc. #149 is in the bubble region, between the large asso- 
ciations #137 and #147, but no nebulosity is seen nearby. 

Assoc. #151 is the association with more OB members of 
our catalogue. It is arc-shaped and part of it follows one bub- 
ble. IC 1613's known supernova remnant (SNR) is included in 
the association. Our catalogue yields 5 detections for the SNR, 
which we discard. We note here one apparently isolated target 
whose position in the Q vs U-B and B-V vs U-B color-color di- 
agrams is the same as for the WO star Several obscured clouds 
fall within the limits of the association. 

The age of the SNR is yet undecided since it shows 
both young and old signatures. While it was origina l ly be- 
lieved to be a bout 2E4 years old (Dod orico & DopitaL Il983t 
iPeimbert et al.LlT98 8) X-ray observations suggest an age of 3E3 



vears dJ ones et al. , 11998 ). The emission lines from the knots re- 



semble those from the Crab Nebula and Cy gnus Loop filaments, 
indic ating that this is a young SNR (Arma ndrofi^ & Masseyl 
11991) . It has been suggested that the SN exploded inside a large 
cavity surrounded by a dense H i cloud. The remnant would have 
just encountered the dense material, pr oducing shocks and hence 
the peculiar features dLozinskava et a l.. 1998, 2003). 

Spectral types are known for three stars in this associa- 
tion. The two bright est blue stars are an 09 I and a Bl la, 
iBresoUn et al.l (l2007l) 's B7 and B6. The authors report strong 
Ha wind emission in the spectrum of B6. The association con- 
tains a red supergi ant, S a ndage & Katem (1976)'s variable star 
[S71]V38 (MO la, iSkifil ^^). The RSG is only 6.7" away 
from association #150. 

Assoc. #154 encloses IC 1613's known WO star The ob- 
ject has very blue Q and U-B colors (V=19.86, Q=-1.35, B- 
V=0.02, U-B=-1.33). The WO is located on the edge of the gi- 
ant H I shell that contains the S W side of the galaxy, and ion- 
izes an asymmetric bipolar H ii region (lAfanas'ev et all l2000t 
[Lozinskava et al., 2001). These authors find at least two extra 
sources in H„ and Hell 4686 narrow-band filters that are com- 
pact knots of gas. Some entries of our OB candidate list are also 
knots of the nebula ionized by the WR. They are characterized 
by very blue Q and red B-V. However, U, B and R VLT-VIMOS 
images suggest that two of the detections (besides the WO) are 
in fact stars. One of them has very blue Q, and the other one, 
very embedded in the WO nebula, has red Q. These stars might 
be members of a cluste r hidden behind the WO, proposed by 
ILoz inskava et al. (2001 ). The b ipolar shell's dynamical age is 
0.3-1E6 years ( Afanas'ev et al.L 12000.) . consistent with the dura- 
tion of the WR stage. 

Assoc. #157: VLT-VIMOS images reveal that the brightest 
star is triple, and unresolved in our catalogue. One member has 
similar colors to the brightest star in #158, which has been found 
to be an LBV-candidate. 



Assoc. #158 is located in the void of the bubble region, 
like #149. One of the members has very red B-V and blue Q. 
Medium resolution spectroscopy soon to be published (Herrero 
et al. 2009, in prep.) has revealed P Cygni profiles of the Balmer 
series and of iron, rendering the target an LBV-candidate. 

Assoc. #161 traces two bubble rims, s haring one of them wi th 
association #147. The brightest target is lBresolin et al.l (l2007l) 's 
All (B9 la). 

Assoc. #162's brightest star is in the intersection with one 
of the bubble rims fo llowed by #161. Its spectral type is 05- 
06 V jBresoUn et all (|2007)'s B2), and it is likely ionizing its 
surrounding nebula. 

Assoc. #175: The association is spatially divided into two 
parts, one at the centre of a bubble, and the other in the intersec- 
tion of such bubble and another that encloses association #185. 
One apparently isolated member has very blue Q ( — 1 .4), and its 
position in the Q vs U-B and B-V vs U-B diagrams is similar to 
the WO. 

Assoc. #176: Two out of its three members are in the centre 
of a small bubble. 

Assoc. #179's brightest star is lBresolin et al.1 ( l2007l) 's A12, a 
B1.5 lab. 

Assoc. #184 is found in the SW outskirts of the galaxy. The 
stars are distributed in a ring shape, around a central red star. 

Assoc. #185: Members are located on the rim of one bubble, 
but they may also be creating the bubble where part of associa- 
tion #175 is located at the SW. The association has a high stellar 
density and all members may expe rience blends. The br ightest 
star is an early B-supergiant (BO la iBresolin et al.l ( 2007l)'s B3: 
B:I: Sandage & KatemI jl976 f)'s B42) although' Lozinskava et al.1 
(2002) classified it as an O supergiant. The study of the latter in 
this region (see comment on association #147 for a description 
of their analysis) found 3 more OB stars: O III, B8 II and one 
with uncertain type 08 III or 04 V . A faint member (V~22) 
with extremely blue Q and U-B, sits in the same position in the 
U-B vs Q diagram as the WO. 

Assoc. #187's brightest member is IBresolin et al.l (l2007l) 's 
A13 a very early O dwarf (03-4 V((f))), which is remarkable 
for a modest association of only 3 members. 

Assoc. #189 is located on the rim of a bubble. The spatial 
distribution suggests that the H ii shell could have been blown 
either by association #185 or by a nearby isolated star 

Assoc. #191 is linear in shape and located in the bubble re- 
gion, close t o association #189. T he brightest star is a late-O 
giant (09 III. lBresoUn eTal] (l2007h 's A15). 

Assoc. #196's brightest star is a Be. iBresoUn et al.1 (l2007l) 's 
A17. Their target A16 (early-B lb) is only 6.1" away. 

7. THE EFFECT OF THE V-FAINT LIMIT ON THE 
FINAL CATALOGUE OF ASSOCIATIONS 

The selection criteria for OB stars were flexible in magnitude 
and colors to admit very extincted stars, since we suspected that 
internal reddening was significant. However, after inspecting the 
associations we have found that the redder faint stars are not 
reddened blue stars, but mostly non-stellar objects or blends. 
We consequently wonder how allowing such a large range of 
V-magnitudes and colors has distorted the resulting catalogue 
of associations. Fainter stars are more abundant, hence closer to 
each other, and may shift the value of Ds to shorter distances. 
By setting Ds to a short value, we may be missing looser groups 
of bright OB candidates. 

We ran the OB association search programme on a subset of 
the initial sample: bright stars with V<21. The resulting list of 
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Fig. 10. Comparison of the results presented in this paper 
(catalogue-A, grey shadows) and the results of the OB associa- 
tion search algorithm run on the same sample of stars but V<21 
(catalogue-B, black). The associations found in both cases agree 
on the bubble region, but the SW side of the galaxy is dominated 
by associations with fainter stars. Catalogue-A misses 3 associ- 
ations of bright stars. 



associations (hereafter catalogue-B) is compared with our results 
(catalogue-A) in Fig.fTOl 

As expected we obtained a longer Ds=H" for the less dense 
group of bright stars. There is good agreement of both cata- 
logues on the bubble NE region of the galaxy. The associations 
in catalogue-B are usually larger and enclose several associa- 
tions of catalogue-A. However, the SW side of the galaxy in the 
H I cavity is apparently dominated by fainter stars, and by limit- 
ing ourselves to the brightest ones we would have missed a large 
number of associations. On the other hand, by enlarging Ds to 
11" we would include three more associations. Whether they are 
real OB associations can only be decided after closer inspection 
of their stellar content. 

We conclude that the criteria used to build catalogue-A 
favoured the identification of OB associations, although some 
extra work was needed afterwards to discard fake identifications. 



8. COMPARISON WITH PREVIOUS CATALOGUES 

Four catalogues of OB associations in IC 1613 have been pub- 
lished to date. The first one was made from visual inspection 
of plates and published by Hodge (1978). More recent works 
byLIvanov ( 1 996h (c a talogu e not public). iGeorgiev et al.l (Il999l) 
and iBorissova et alj (|2004|) have used an automatic procedure 
very similar to ours. Our work has produced the most numerous 
catalogue: a total of 163 stellar groups, almost triplicating the 
60 associations fou nd by Borisso va et al. Hodge ( 197 8) listed 
20 ass ociations and llvanovl (Il996h only found 6. lOeorgiev et alj 
(119991) studied a smaller field. 



Fig. 11. Comparison of the association s fou nd in this work 
(grey), with those of iBorissova et al.l ( 120041) (black). The 
quadrants mark the same positions as in Fig. |9] The over- 
all distribution of the associations in both catalogues agrees. 
Considering the associations in the H i cavity, our results sep- 
arate Borissova et al. (2004) group s at the SB into many pieces. 
At the NW end, IBorissova et al.l (l2004l) find 4 associations 
(G51,G56,G57,G58) that do not have a counterpart in our cat- 
alogue. Other unmatched associations are their G26 and G55. 



Direct co mparison of our results with the seminal work of 
lHodgel(ll978h can only be done qualitatively, since his catalogue 
is marked on low resolution photographic plates of the galaxy. 
The associatio ns we have found follow the same distribution as 
'Hodgd (Il978| )'s, but our c ode breaks them into many smaller 
ones. Georgie v et al.l ( Il999h observed the same effect in their re- 
sults for the bubble area of the galaxy, and pointed out that their 
associations are the bright cores of Hodge's. We also find new 
associations both at the core of IC 1613 and the periphery. Both 
facts are explained by the higher sensitivity and finer spatial res- 
olution of this work. 

The associations found by Georgiev et al.l (11999 *) are pro- 
vided on detector coordinates or plotted over a poor qual- 
ity image, and since the same team studied a wider field in 
IBorissova et alj ( |2004|) we preferred the latter for comparison. 
We use the coordinates and radius of the associations in their 
Table 1 to compare with ours (see Fig.irTTl. We noted a discrep- 
ancy between the centres provided in their list and the position 
of some associations they marked over a picture of IC 1613 in 
their Fig. 1 . By providing electronically the list of OB members 
of each association, we expect to minimize these kind of incon- 
sistencies and to make direct inter-comparisons easier. 

Th e overall distribution of our catalogue and lBorissova et akl 
(l2004l) 's associations agrees. The data they used have a plate 
scale similar to ours but again, our catalogue separates most of 
their associations into two or more, especially in the SW part of 
the galaxy enclosed in the H i cavity (see Fig.lTTli. Our catalogue 
has additional associations in the outerreaches of IC 1613, as ex- 
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pected because of our extended field of view, but also in the over- 
lapping areas. The likely reason why these associations were not 
registered by Borissova et al. is that they do not reach the min- 
imum number of elements considered in their work (Nmin=4). 
In other cases the members are faint, beyond Borissova et al.'s 
limiting V-magnitude. These authors found 3 associations (G26, 
G51, G55) far from any of ours (see Fig. [TTI ): we do find OB 
stars on those locations, but they are separated more than Ds and 
hence not included in our catalogue. There is also a mismatch be- 
tween associations G57-G58 and the centre of our counterparts. 

Th e different number of associations found by Hodge 
(Il978h . [Borissova et al. ' (2004) and ourselves may be iDctter ex- 
plained in light of the experiment described in Sect. Q Hodge 
only considered stars with V<21 to find associations, and 
Borissova et al. have completeness problems for fainter magni- 
tudes. Catalogue-B consists of 31 associations o f average diame- 
ter s imilar to the mean size of the associations of lBorissova et alj 
(12004.) . decreasin g the discrepanc i es. Ho wever, catalogue-B still 
does not include iBorissova et al.l (l2004l) 's associations without 



counterpart in catalogue-A. 



9. SUMMARY AND FUTURE WORK 

This paper presents a deep multi-band photometric catalogue of 
IC 1613 with high astro me trie accuracy. Using the uni vocal re- 
lation between the reddening-free Q parameter and OB spectral 
types, we have built a list of candidate young massive stars. With 
a friends-of-friends algorithm developed by our group, we have 
studied how they cluster into associations. A total of 163 associ- 
ations were found (after discarding 5 galaxies and 29 cases with 
less than Nmin valid members), ranging from small groups of 
3 stars to large complexes of more than 20 members in the cen- 
tre of H II shells. 24 associations have 10 or more OB candidates. 
Our findings are dependent on the choice of the two free parame- 
ters of the friends-of-friends algorithm: the search distance (6") 
and the minimum number of OB stars defining an association 
(3). 

The resulting set of candidate OB stars is ready to be ob- 
served with multi-object spectrographs. To achieve the required 
astrometric precision, a new solution for the WFC at the INT 
was developed. 

The photometric analysis leaves an interesting by-product, 
an extinction map of IC 1613. Most of the OB stars exhibit a 
color excess larger than the foreground value. Extinction is not 
uniform and does not have a smooth distribution. The points of 
enhanced reddening are spread throughout the galaxy and show 
a slightly higher concentration at the edge of the H i cavity that 
contains the SW side of the galaxy. 

The positive confrontation of our list of candidate OB stars 
with published spectral types, endorse our target selection crite- 
ria based mainly on Q and secondarily on U-B colors. Four stars 
are found in the same position of the U-B vs Q diagram as the 
known Wolf-Rayet star of the galaxy. 

This article is the first effort towards a comprehensive de- 
scription of the massive stellar population of IC 1613, as a proxy 
for the low-metallicity early-Universe. The next step is the char- 
acterization of the physical properties of the associations pre- 
sented here, especially mass and age. The study of the color- 
magnitude and color-color diagrams of the catalogue associa- 
tions and their comparison with theoretical isochrones will be 
presented in part-II. It will provide interesting clues about star 
formation processes of massive stars and their IMF in low metal- 
licity environments. 
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Appendix A: CATALOGUE OF OB ASSOCIATIONS 

A.1. Complete Catalogue 

The output list of OB associations provided by our code is pre- 
sented in Table lA. 13] 

The content of different columns is: (1) OB association 
identification number. (2) and (3) Right Ascension (RA) and 
Declination (DEC) of the centre determined from the average of 
the positions of the OB members, equinox J2000.0. (4) Number 
of OB members (Nqb) of the association. (5) Total number of 
members (Ntot); this figure may not be representative since our 
photometric selection criteria is biased against red objects by de- 
manding high U-magnitude accuracy. (6) Quadrant allocation in 
Fig-El (7) Notes: 'Gal' - the association is actually a galaxy; 'd' 
- association discarded because it does not have enough valid 
members (i.e. visual inspection reveals that some are blends or 
gas knots); '?' - dubious association (associations with a small 
number of members, with a significant fraction under suspicion 
of being blends or not blue stars); 'n' - an entry in Sect. 16. II 

A.2. Member Listing Example: Association #16 

Our code provides us with a list of the OB candidate members 
of each association and also the red population. Defining OB as- 
sociations by their members, rather than by hand-drawn lines on 
photographic plates, allows a direct comparison of results be- 
tween different works and with observations. The member list- 
ings of all the associations of our catalogue will be available 
on-line. 

As an example, we list the members of association #16 in 
Table IATT41 The content of different columns is: (1) Star iden- 
tification number in our catalogue. (2) and (3) Right Ascension 
and DecUnation, equinox J2000.0. (4), (6), (8), (10) and (12) 
magnitudes. (5), (7), (9), (11) and (13) error of magnitude mea- 
surements. (14) and (15) DAOPHOT's;^' and # parameters. (16) 
Identification number in the 2MASS catalogue. 

A.3. Finding Charts of the Associations 

The positions of the associations were marked on the galaxy in 
Fig. |9] Zooms into different galactic sectors are provided in this 
Appendix, Figs. lA.HIA.9l 



M. Garcia et al.: OB Associations in IC 1613, Online Material p 2 
Table A.13. Catalogue of OB associations in IC 1613. 



ID 


RA [deg] 
(J2000.0) 


DEC [deg] 
(J2000.0) 


Nob 




Quad 


Notes 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


1 


16.096912 


2.052357 


8 


10 


3 


Gal 


2 


16.102215 


2.110539 


3 


3 


2 




3 


16.103825 


2.152657 


4 


4 


1 




4 


16.103556 


2.127854 


3 


3 


2 


d 


5 


16.108092 


2.114839 


3 


3 


2 


? 


6 


16.110625 


2.110028 


3 


3 


2 




7 


16.112182 


2.139603 


4 


4 


2 




8 


16.116793 


2.148229 


4 


4 


1,2 




9 


16.116749 


2.181977 


5 


5 


1 




10 


16.120176 


2.176666 


3 


3 


1 


n 


11 


16.125909 


2.149869 


4 


4 


1 




12 


16.125697 


2.104176 


3 


3 


2 


d 


13 


16.127069 


2.128816 


4 


4 


2 




14 


16.129429 


2.097499 


4 


4 


2 




15 


16.130444 


2.093961 


3 


3 


2 




16 


16.133370 


2.165873 


10 


12 


1 


n 


17 


16.130397 


2.059373 


3 


3 


3 




18 


16.132333 


2.150075 


4 


4 


1,2 




19 


16.132298 


2.169800 


3 


4 


1 


? 


20 


16.133645 


2.077595 


4 


4 


3 




21 


16.135528 


2.161271 


17 


18 


1 


n 


22 


16.138624 


2.161875 


3 


3 


1 




23 


16.140158 


2.098123 


4 


4 


2 


7 


24 


16.143162 


2.105718 


5 


5 


2 




25 


16.150666 


2.098182 


7 


9 


2 


d 


26 


16.149934 


2.104582 


3 


3 


2 




27 


16.150782 


2.094274 


5 


5 


2 


n 


28 


16.154294 


2.093873 


4 


4 


2 


n 


29 


16.156629 


2.107164 


3 


3 


2 


? 


30 


16.157865 


2.118325 


4 


4 


2 




31 


16.160063 


2.086166 


3 


4 


2 




32 


16.165508 


2.129537 


7 


7 


2 


Gal 


33 


16.167062 


2.145645 


5 


7 


2,5 




34 


16.173860 


2.091680 


7 


8 


5 




35 


16.172452 


2.151775 


8 


9 


4 




36 


16.173835 


2.146373 


4 


7 


4,5 




37 


16.174840 


2.174687 


4 


5 


4 




38 


16.174955 


2.083569 


3 


3 


5 




39 


16.176810 


2.138951 


3 


3 


5 


d 


40 


16.177356 


2.153791 


3 


3 


4 


? 


41 


16.181246 


2.070946 


4 


4 


6 




42 


16.180837 


2.156761 


3 


3 


4 


d 


43 


16.182817 


2.112977 


3 


5 


5 


n 


44 


16.186156 


2.108783 


16 


17 


5 




45 


16.185969 


2.081011 


3 


3 


5,6 


? 


46 


16.187206 


2.165141 


3 


3 


4 


d 


47 


16.187683 


2.158649 


3 


5 


4 


? 


48 


16.190017 


2.110850 


9 


9 


5 


n 


49 


16.189936 


2.130523 


3 


5 


5 


d 


50 


16.190243 


2.163109 


3 


5 


4 


d 


51 


16.194281 


2.106449 


8 


10 


5 


n 


52 


16.193448 


2.127703 


3 


4 


5 


d 


53 


16.194693 


2.132499 


7 


8 


5 


n 


54 


16.198802 


2.111230 


9 


10 


5 


n 


55 


16.197736 


2.101776 


3 


3 


5 


d 


56 


16.203010 


2.107146 


16 


19 


5 




57 


16.200013 


2.100085 


4 


4 


5 


? 


58 


16.200024 


2.125184 


6 


6 


5 


? 


59 


16.200528 


2.133959 


3 


4 


5 


d 


60 


16.202105 


2.078538 


6 


6 


5,6 


? 


61 


16.202044 


2.128357 


4 


4 


5 


n 


62 


16.202794 


2.102484 


3 


3 


5 


d 


63 


16.204786 


2.112106 


7 


7 


5 




64 


16.204397 


2.124905 


5 


5 


5 


_ 


65 


16.204474 


2.095443 


4 


5 


5 


d 


66 


16.205571 


2.130890 


8 


10 


5 


? 


67 


16.206485 


2.126462 


4 


5 


5 


? 


68 


16.206628 


2.099115 


3 


3 


5 


? 


69 


16.209295 


2.117948 


6 


8 


5 


n 


70 


16.208526 


2.094393 


3 


3 


5 


d 


71 


16.208783 


2.083076 


3 


4 


5 


d 


72 


16.209589 


2.101978 


3 


3 


5 


7 
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Table A.13. continued. 



ID 


RA [deg] 
(J2000.0) 


DEC [deg] 
(J2000.0) 


Nob 


Ntot 


Quad 


Notes 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


73 


16.209537 


2.131091 


5 


5 


5 


- 


74 


16.210802 


2.111945 


3 


3 


5 


d 


75 


16.211725 


2.106787 


5 


8 


5 


?n 


76 


16.213125 


2.076123 


12 


13 


6 


- 


77 


16.211226 


2.072083 


4 


4 


6 


? 


78 


16.212151 


2.158723 


4 


4 


4 


? 


79 


16.211941 


2.127620 


3 


3 


5 


? 


80 


16.212032 


2.104437 


3 


5 


5 


_ 


81 


16.213312 


2.131232 


3 


3 


5 


- 


82 


16.213018 


2.124339 


4 


5 


5 


? 


83 


16.214367 


2.061752 


3 


3 


6 


n 


84 


16.216014 


2.101695 


10 


10 


5 


- 


85 


16.213772 


2.118912 


3 


3 


5 


?n 


86 


16.214380 


2.126652 


3 


3 


5 


?n 


87 


16.215442 


2.115048 


5 


5 


5 


- 


88 


16.217727 


2.107367 


16 


17 


5 


- 


89 


16.218450 


2.183555 


6 


6 


4 


- 


90 


16.218416 


2.163930 


3 


3 


4 


? 


91 


16.219528 


2.079689 


6 


6 


5,6 


- 


92 


16.219961 


2.098888 


3 


3 


5 


d 


93 


16.220890 


2.084779 


6 


7 


5 


- 


94 


16.219702 


2.103431 


4 


4 


5 


? 


95 


16.220188 


2.149358 


3 


3 


4 


- 


96 


16.222186 


2.096009 


4 


5 


5 


n 


97 


16.222757 


2.103359 


3 


3 


5 


n 


98 


16.223875 


2.146328 


6 


7 


4,5 


n 


99 


16.223379 


2.058472 


4 


4 


6 


n 


100 


16.224860 


2.169550 


6 


6 


4 


- 


101 


16.223822 


2.179711 


3 


3 


4 


d 


102 


16.225417 


2.089271 


3 


3 


5 


d 


103 


16.229149 


2.096880 


19 


23 


5 


n 


104 


16.226815 


2.155823 


7 


10 


4 


?n 


105 


16.226303 


2.137465 


3 


3 


5 


dn 


106 


16.227180 


2.101775 


14 


18 


5 


n 


107 


16.227032 


2.107051 


3 


3 


5 


d 


108 


16.227880 


2.084607 


4 


6 


5 


? 


109 


16.227894 


2.147403 


5 


5 


4,5 


? 


110 


16.228626 


2.173969 


3 


4 


4 


? 


111 


16.229030 


2.104253 


5 


5 


5 


? 


112 


16.231603 


2.093674 


3 


4 


5 


? 


113 


16.231901 


2.103999 


5 


6 


5 


- 


114 


16.235155 


2.082594 


22 


23 


5,6 


n 


115 


16.233123 


2.089860 


3 


3 


5 


d 


116 


16.236758 


2.089753 


6 


7 


5 


n 


117 


16.237337 


2.079615 


3 


4 


5,6 


n 


118 


16.238855 


2.135455 


3 


3 


5 


d 


119 


16.238087 


2.173274 


3 


3 


4 


- 


120 


16.240415 


2.096584 


7 


7 


5 


n 


121 


16.240945 


2.177397 


4 


4 


4 


- 


122 


16.242760 


2.079217 


15 


19 


5,6,8,9 


n 


123 


16.241480 


2.089006 


14 


15 


5,8 


- 


124 


16.241957 


2.084267 


7 


8 


5,8 


n 


125 


16.242312 


2.124345 


5 


5 


5,8 


n 


126 


16.243098 


2.093806 


5 


5 


5,8 


n 


127 


16.248152 


2.153727 


49 


54 


4,7 


n 


128 


16.242096 


2.149247 


4 


5 


4,5 


n 


129 


16.242767 


2.101736 


6 


6 


5,8 


- 


130 


16.244092 


2.158721 


4 


5 


4,7 


_ 


131 


16.243844 


2.098030 


3 


3 


5,8 


- 


132 


16.248173 


2.090002 


13 


14 


8 


- 


133 


16.247413 


2.085148 


5 


5 


8 


- 


134 


16.246209 


2.071744 


3 


3 


9 


n 


135 


16.249324 


2.097876 


31 


37 


8 


- 


136 


16.247844 


2.148337 


4 


4 


7,8 


n 


137 


16.253227 


2.178541 


46 


48 


7 


n 


138 


16.248630 


2.163478 


6 


6 


7 


? 


139 


16.248756 


2.136402 


3 


3 


8 


d 


140 


16.248518 


2.144179 


3 


3 


8 


?n 


141 


16.249494 


2.081922 


3 


3 


8 


Gal 


142 


16.250423 


2.054868 


4 


4 


9 




143 


16.250916 


2.170979 


3 


3 


7 




144 


16.251720 


2.146901 


3 


3 


8 


n 


145 


16.252874 


2.086635 


9 


9 


8 
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Table A.13. continued. 



ID RA[deg] DEC[deg] Nqb Ntut Quad Notes 
(J2000.0) (J20{)0.0) 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


146 


16.254549 


2.104795 


10 


10 


8 


n 


147 


16.257084 


2.158485 


53 


58 


7 


n 


148 


16.254275 


2.092485 


6 


7 


8 


15 


149 


16.254539 


2.171786 


3 


3 


7 


?n 


150 


16.255150 


2.142742 


5 


6 


8 




151 


16.260516 


2.140498 


99 


103 


7,8 


n 


152 


16.256125 


2.190009 


4 


4 


7 




153 


16.256249 


2.125097 


5 


5 


8 




154 


16.256969 


2.072318 


11 


12 


9 


n 


155 


16.258348 


2.100803 


9 


10 


8 




156 


16.258377 


2.096784 


3 


5 


8 


d 


157 


16.259700 


2.153904 


4 


4 


7 


?n 


158 


16.258826 


2.173986 


3 


3 


7 


?n 


159 


16.258869 


2.029636 


3 


3 


9 


Gal 


160 


16.259788 


2.106187 


5 


5 


8 




161 


16.262861 


2.161338 


11 


12 


7 


n 


162 


16.262323 


2.169053 


3 


5 


7 


n 


163 


16.262203 


2.084492 


3 


3 


8 




164 


16.262549 


2.196617 


3 


4 


7 


? 


165 


16.264022 


2.181578 


4 


4 


7 




166 


16.264283 


2.050493 


6 


6 


9 


Gal 


167 


16.266025 


2.101659 


3 


3 


8 


? 


168 


16.265918 


2.097347 


3 


4 


8 


? 


169 


16.266384 


2.192655 


4 


4 


7 




170 


16.266803 


2.166934 


3 


3 


7 


? 


171 


16.266873 


2.138073 


4 


5 


8 




172 


16.267850 


2.178595 


3 


3 


7 


d 


173 


16.268078 


2.152881 


3 


3 


7 




174 


16.268225 


2.203677 


3 


3 


7 




175 


16.270523 


2.157452 


28 


30 


7 


n 


176 


16.268470 


2.163499 


3 


3 


7 


n 


177 


16.271771 


2.186759 


5 


5 


7 




178 


16.271105 


2.200210 


4 


6 


7 




179 


16.271636 


2.142294 


3 


3 


8 


?n 


180 


16.273082 


2.147837 


4 


7 


7,8 




181 


16.272820 


2.084834 


4 


4 


8 


d 


182 


16.273732 


2.204206 


4 


4 


7 




183 


16.274270 


2.150566 


3 


4 


7 


d 


184 


16.273774 


2.177004 


3 


3 


7 


?n 


185 


16.276164 


2.158762 


17 


17 


7 


n 


186 


16.275929 


2.194830 


3 


3 


7 


? 


187 


16.276569 


2.178761 


3 


3 


7 


n 


188 


16.277759 


2.145682 


4 


5 


8 


? 


189 


16.278765 


2.164795 


4 


4 


7 


?n 


190 


16.280133 


2.183874 


3 


3 


7 


d 


191 


16.283802 


2.166530 


13 


13 


7 


n 


192 


16.286930 


2.151835 


3 


3 


7 


? 


193 


16.288172 


2.181669 


4 


4 


7 


? 


194 


16.290408 


2.170542 


6 


7 


7 


? 


195 


16.290153 


2.185749 


4 


4 


7 


? 


196 


16.295735 


2.207267 


5 


5 


7 


n 


197 


16.306976 


2.139009 


8 


8 


8 




198 


16.316533 


2.200474 


3 


3 


7 


d 



Table A.14. Members of association #16. 



OB-stars 


in Association 




























ID 


KA [degj 


UnL. [clegj 


V 


cr-V 


B 


cr-B 


I 


cr-I 


R 


C7--R 


u 


C7--U 


y 

A. 


# 


ID.2MASS 




(J ZUUU.U) 


(J ZUUU.U; 




























(1) 






(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


17183 


lo.l299o4 


2.165500 


19.885 


0.0074 


19.827 


0.0055 


19.894 


0.0143 


19.936 


0.0089 


19.220 


0.0121 


1.405 


0.051 




17397 




Z.loo5o3 


22.656 


0.0477 


22.636 


0.0345 


23.124 


0.1523 


22.808 


0.0642 


21.689 


0.0398 


1.219 


0.531 




17512 


16.131291 


2.166572 


21.861 


0.0257 


21.716 


0.0182 


22.190 


0.0629 


21.964 


0.0351 


20.852 


0.0195 


1.244 


0.322 




17896 


16.132740 


2.166766 


21.819 


0.0222 


21.655 


0.0162 


21.958 


0.0546 


21.987 


0.0335 


20.550 


0.0181 


1.200 


0.271 




17945 


16.132886 


2.166397 


22.193 


0.0325 


22.129 


0.0245 


21.988 


0.0672 


22.139 


0.0372 


21.379 


0.0335 


1.259 


0.429 




18051 


16.133337 


2.165283 


22.313 


0.0301 


22.172 


0.0209 


22.379 


0.0962 


22.485 


0.0511 


21.324 


0.0289 


1.193 


0.380 




18159 


16.133736 


2.164780 


21.968 


0.0357 


21.943 


0.0200 


21.942 


0.0957 


22.053 


0.0567 


21.095 


0.0286 


1.629 


0.776 




18444 


16.134765 


2.165754 


22.547 


0.0351 


22.455 


0.0254 


22.586 


0.1204 


22.652 


0.0590 


21.678 


0.0350 


1.158 


0.472 




18867 


16.136389 


2.165391 


20.818 


0.0132 


20.641 


0.0090 


21.012 


0.0270 


20.981 


0.0164 


19.598 


0.0139 


1.427 


0.226 




19245 


16.137705 


2.165704 


22.499 


0.0302 


22.462 


0.0236 


22.531 


0.0812 


22.493 


0.0416 


21.912 


0.0416 


0.999 


0.037 





Non OB-stars in Association 



ID 


RA [deg] 


DEC [deg] 


V 


0--V 


B 


cr-B 


I 


0--I 


R 


0--R 


U 


0--U 


X 


# 


ID.2MASS 




(J2000.0) 


(J2000.0) 




























(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


16825 


16.128596 


2.164702 


21.332 


0.0173 


21.617 


0.0150 


20.855 


0.0280 


21.120 


0.0200 


21.678 


0.0367 


1.266 


0.383 




18116 


16.133573 


2.166236 


21.772 


0.0352 


22.151 


0.0217 


21.064 


0.0657 


21.464 


0.0485 


22.176 


0.0466 


1.789 


0.661 
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1h4m22s 



1h4m28s - 



1h4m34s - 



1 h4m40s 



1. 

• 

■ 

« 


1 


♦ 

a 

- • • <<b 

i . * ^ 


il r 

• 






1- 

• 


■ . ■ ■ -V 


1 


-^22 ' 

• 


* 

« 

■ ■ 


■ 




• 

• 

* 


* ' • 

1 • . . 





2d12m59s 



2d1 1m37s 



2d 10ml 5s 



2dSm53s 



Fig. A.l. INT-WFC U-band image of a portion of IC 1613 corresponding to quadrant 1 of Fig. |9] The OB associations located in 
this quadrant are indicated by a black contour. Their numbers are: #3, #8-#l 1, #16, #18, #19, #21 and #22. 
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1h4m22s 



1h4m28s 



1h4m34s - 



1 h4m40s 




2d8m5Js 



2d7m32s 



2d6m10s 



2d4m48s 



Fig. A.2. Same as\AJ] quadrant 2. Associations present in this quadrant: #2, #4-#8, #12-#15, #18, #23-#33. 
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1h4m22s 



1h4m2as - 



1 h4m34s - 



1 h4m405 







1 






« * 


m 


■ 

• 


* - 

♦ 

■ 


m 




• 

• 


« 

• 

• 


M m 


• 

. 3 

• 


* 

« 


• 

• • 

• 

* 

• 

• • 

1 



2d4m4Ss 2d3m26s 2d2m4g 

Fig. A.3. Same as lA.ll quadrant 3. The enclosed associations are #1, #17 and #20. 



2d0m425 
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1 h4m40s 



1 h4m46s ~ 



1 h4m59s 



* 4 



2d12m59g 



1. 



36 



l37 




35 , 
Q42 



■ '^9 



0,110 ^ 



0^121 




95/3 
98 

104 

109 



.126 



2d1 1m37s 



2d 1 0ml 5s 



2d8m53s 



Fig. A.4. Same as|Aj] quadrant 4. Associations present in this quadrant: #35-#37, #40, #42, #46, #47, #50, #78, #89, #90, #95, #98, 
#100, #101, #104,# 109#, 1 10, #1 19, #121, #127, #128 and #130. 
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1 h4m40s 



1 h4m46s - 



1h4m52s - 



1 h4m59s 




2d8m55s 



2d7m31s 



2d6m10s 



2d4m4S5 



Fig.A.5. Same as|Aj] quadrant 5. Associations present in this quadrant: #33, #34, #36, #38, #39, #43-#45, #48, #49, #51-#75, 
#79-#82, #84-#88, #91-#94, #96-#98, #102-#103, #105-#109, #1 1 1-#1 18, #120, #122-#126, #128, #129 and #131. 
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1 h4m40s 



1 h4m46s - 



1h4m52s 



1 h4m59s 




2d4m4Ss 



2d3m26s 



2d2m4s 



2d0m425 



Fig, A.6. Same as\AA\ quadrant 6. The enclosed associations are #41, #45, #60, #76, #77, #83, #91, #99, #1 14, #1 17 and #122. 
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Fig.A.7. Same as [All quadrant 7. Associations present in this quadrant: #127, #130, #136-#138, #143, #147, #149, #151, #152, 
#157, #158, #161, #162, #164, #165, #169, #170, #172-#178, #180, #182-#187, #189-#196 and #198. 
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Fig.A.8. Same as|A3] quadrant 8. Associations present in this quadrant: #122-#126, #129, #131-#133, #135, #136, #139-#141, 
#144-#146, #148, #150, #151, #153, #155, #156,# 160, #163, #167, #168, #171, #179-#181, #188 and #197. 
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1h4m59s 



1h5m5s 



IhSmlls - 



1h5m17s 




2d4m48s 



2d3m26s 



2d2m4s 



2d0m423 



Fig, A.9. Same as\AA\ quadrant 9. The enclosed associations are #122, #134, #142, #154,# 159 and #166. 



